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ABSTRACT 

We outline a scenario of galaxy formation in which the gas in galaxy-forming regions 
was preheated to high entropy by vigorous energy feedback associated with the for¬ 
mation of stars in old ellipticals and bulges and with AGN activity. Such preheating 
likely occurred at redshifts z 2 - 3, and can produce the entropy excess observed 
today in low-mass clusters of galaxies without destroying the bulk of the Lya forest. 
Subsequent galaxy formation is affected by the preheating, because the gas no longer 
follows the dark matter on galaxy scales. The hot gas around galaxy haloes has very 
shallow profiles and emits only weakly in the X-ray. Cooling in a preheated halo is not 
inside-out, because the cooling efficiency does not change significantly with radius. 
Only part of the gas in a protogalaxy region can cool and be accreted into the final 
galaxy halo by the present time. The accreted gas is likely in diffuse clouds and so 
does not lose angular momentum to the dark matter. Cluster ellipticals are produced 
by mergers of stellar systems formed prior to the preheating, while large galaxy disks 
form in low-density environments where gas accretion can continue to the present 
time. 

Key -words: galaxies: formation - galaxies: structure - galaxies: spiral - galaxies: 
elliptical - galaxies: clusters 


1 INTRODUCTION 

In the past 20 years, the Cold Dark Matter (CDM) cos¬ 
mogony (Peebles 1982; Blumenthal et al. 1984) has become 
the most appealing scenario for the formation of structure 
in the universe. In the CDM cosmogony, the universe is as¬ 
sumed to be dominated by CDM, and the structure in the 
universe forms in a hierarchical fashion, with dark matter 
particles aggregating into larger and larger clumps (dark 
haloes) in the passage of time. In this scenario, galaxies are 
assumed to form as gas cools and condenses in dark haloes 
(White & Rees 1978), and the properties of galaxies are de¬ 
termined by the properties of their host haloes, such as mass, 
density profile, angular momentum, and merging history. 
This standard scenario of galaxy formation has been quite 
successful in predicting many of the observed properties of 
the galaxy population, such as disk sizes and kinematics 
(e.g. Fall & Efstathiou 1980; Dalcanton, Spergel & Summers 
1997; Mo, Mao & White 1998, hereafter MMW; Avila-Reese, 
Firmani & Hernandez 1998; Heavens & Jimenez 1999; Mo & 
Mao 2000; Navarro & Steinmetz 2000; van den Bosch 2000), 
and the distributions of luminosity, colour and morphology 
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(e.g. White & Frenk 1991; Kauffmann, White & Guiderdoni 
1993; Cole et al. 1994; Somerville & Primack 1999). 

One key ingredient in the standard scenario is a mech¬ 
anism which can prevent gas from cooling too fast in dark 
matter haloes, because without such a mechanism, we may 
run into the following set of problems. 

• The overcooling problem: In a hierarchical model, all 
dark matter was in small haloes at high redshifts. Since ra¬ 
diative cooling of gas in these haloes is effective, all the gas 
would have cooled and formed stars, leaving no gas for pro¬ 
ducing quasar absorption line systems and for the forma¬ 
tion of galaxy disks at low redshifts. This is the ‘overcooling 
problem’ first pointed out by White & Rees (1978). 

• The baryon-fraction problem: Based on cosmic nucle¬ 
osynthesis (e.g. O’Meara et al. 2001) and Cosmic Microwave 
Background Radiation (e.g. de Bernardis et al. 2002), the 
baryon density in the universe is Db,o ~ 0.02h~^. For a uni¬ 
verse with flo ~ 0.3 and h « 0.7, the implied overall baryon 
fraction is about /b = flB,o/flo ~ 0.13, consistent with the 
baryon fraction in rich clusters (e.g. White et al. 1993; Ettori 
& Fabian 1999). On the other hand, modelling of galaxy disk 
formation in CDM haloes showed that the gas mass which 
settles into the disk is only ~ 5% of the total halo mass and 
a much higher fraction would produce disks which might be 
unstable and have too peaked rotation curves (e.g. MMW). 
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Since the total amount of hot gas in a galaxy halo is small, 
a consistent model of disk formation requires that only part 
of the gas in a protogalaxy region ends up in the final halo. 
The question is then what regulates the baryon fraction in 
galaxy haloes. 

• Angular-momentum problem: Since gas can cool 
rapidly in small progenitor haloes to form condensed cen¬ 
tral objects, it can lose most of its angular momentum to 
the dark matter during galaxy assembly (e.g. Navarro & 
White 1994). This leads to serious problems when compar¬ 
ing with real galaxies: the disks formed in this way are too 
small. 

• X-ray halo problem: In the standard scenario, gas in a 
galactic halo is assumed to be shock-heated to the virial tem¬ 
perature (which is about 10® K) by gravitational collapse, 
and (massive) galactic haloes are predicted to emit in X- 
ray, with a bolometric luminosity of the order 10“^^ ergs“^. 
The predicted luminosity is about one order of magnitude 
higher than what is observed in the haloes of spiral galax¬ 
ies (e.g. Benson et al. 2000). This is a well-known problem 
(e.g. White & Frenk 1991) but has been largely ignored in 
theoretical studies of galaxy formation so far. 

Two mechanisms have been proposed to solve these 
problems. The first is photoionization heating of the in- 
terglactic medium (IGM) by the UV background. Detailed 
analyses showed that photoionization heating is effective in 
prohibiting gas cooling only in small haloes, with circular ve¬ 
locity Vc ^ 50kms“^ (e.g. Efstathiou 1992; Gnedin 2000). 
Thus, while photoionization heating may be able to keep 
some of the IGM in a diffuse form (and so alleviating the 
overcooling problem), its effect on gas assembling and cool¬ 
ing in galactic haloes may not be sufficient to solve the other 
problems mentioned. The second possibility is heating by the 
energy feedback from star formation (Larson 1974; White & 
Rees 1978; Dekel & Silk 1986). In many semi-analytic models 
of galaxy formation (e.g. White & Frenk 1991; Kauffmann et 
al. 1993; Cole et al. 1994; Somerville & Primack 1999; Kauff¬ 
mann et al. 1999), this feedback is implemented in such a 
way that the gas in galactic haloes is assumed to be effec¬ 
tively heated by supernova explosions associated with the 
star formation in the galaxy, and the heated gas is assumed 
either to be retained within the halo or ejected from the halo 
to some larger scale. Detailed modelling showed that the to¬ 
tal energy feedback from supernova explosions in a galaxy is 
sufficient to keep a large amount of gas at a high tempera¬ 
ture, and the feedback efficiency can be tuned to reproduce 
many of the observed properties of the galaxy population. 
Unfortunately, our understanding of such feedback is still 
very limited, and it is unclear how it operates in detail. 

In this paper, we consider another possible implemen¬ 
tation of the feedback process. We assume that the IGM in 
galaxy-forming regions was heated to some high entropy at 
some high redshift during an early period of vigorous star 
formation and associated AGN activity. We examine the ef¬ 
fect of such preheating on subsequent galaxy formation. This 
implementation is motivated by three lines of observational 
evidence, (i) Observations show that many stars in the uni¬ 
verse may have formed at z ^ 2 in systems reminiscent of 
local starburst galaxies (see Heckman 2001 and references 
therein). Such galaxies can produce supernova-driven out¬ 
flows which may be able to heat the IGM in the surrounding 


regions, (ii) The comoving number density of AGNs is high¬ 
est around redshifts 2 = 2-3 (Shaver et al. 1996). Since a lot 
of kinetic energy can be released from AGNs, the surround¬ 
ing IGM can be heated to a high temperature, if some of this 
kinetic energy is thermalized (e.g. Inoue & Sasaki 2001). (iii) 
Observations of X-ray clusters and groups suggest that there 
is an entropy excess in these systems relative to that ex¬ 
pected from gravitational collapse alone (e.g. Ponman, Gan¬ 
non & Navarro 1999; Lloyd-Davies et al. 2000). This can 
be explained if the IGM in protocluster (protogroup) re¬ 
gions was preheated to a high temperature (e.g. Ponman et 
al. 1999; Tozzi & Norman 2001). The implementation con¬ 
sidered here is related to some of the early implementations 
(e.g. Nulsen & Fabian 1995, 1997; Sommer-Larsen, Gelato & 
Vedel 1999), but has some distinct properties. We emphasize 
the importance of early starbursts and associated AGN ac¬ 
tivity, and we consider heating of the gas while it was still in 
diffuse form outside galactic haloes. Furthermore, our model 
presents more detailed discussions about the properties of 
the preheated gas in connection to the puzzling features of 
observed galaxies mentioned above. 

Preheating affects subsequent galaxy formation because 
the heated gas cannot collapse into small dark haloes and 
so no longer follows the bottom-up clustering hierarchy of 
the dark-matter component. We show that the introduc¬ 
tion of preheating can help to solve some of the problems 
mentioned above. The paper is structured as follows. In Sec¬ 
tion we describe our assumptions on the properties of the 
preheated IGM and analyze gas cooling and accretion onto 
dark haloes in such a medium. The formation of galaxies in 
the present scenario is described in Section In Section ^ 
we discuss briefly how the IGM may have been preheated. 
Further discussion of the present model and a summary of 
our main conclusions are given in Section Throughout 
this paper, when necessary, we adopt the standard ACDM 
cosmogony, with matter density flo = 0.3, a cosmological 
constant corresponding to Da = 0.7, and a normalization 
(78 = 0.9 for the power spectrum. We write the Hubble con¬ 
stant as Ho = 100hkms“^Mpc“^ and take h = 0.7. 

2 GAS ACCRETION BY DARK HALOES IN A 

PREHEATED INTERGALACTIC MEDIUM 

2.1 The preheated IGM 

For a (completely ionized) cosmic gas, the mean number 
density of electrons at redshift 2 is 

rie « 2.0 X 10"’^ (1-b 2)® cm“® , (1) 

where Db,o is the cosmic density parameter of baryons, 
and ~ 0.02 as given by cosmic nucleosynthesis (e.g. 

O’Meara et al. 2001). Now consider a volume in which the 
electron density is rie = ne{l + d) and the temperature is T. 
We define an entropy parameter as 


where it is conventional to express the temperature in units 
of keV (IkeV corresponds to T « 1.2 x 10^ K). If the gas is 
adiabatic, then S is conserved and the temperature of the 
gas can be written in terms of S as 
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r = 4.4 X 10'‘<Sioo(l + + zfK, 


(3) 


where 


<5ioo = 


( ^ \ 


Vl00h-i/3keVcm2/ 

1 0.02 ) 


2/3 


- 1/3 


(4) 


ho .7 = hjO.r. Note that the h dependence at the end of the 
equation appears because the observed value of S has such a 
dependence. For a given S the required temperature is lower 
if the gas density is lower. For example, if <Sioo = 1 and if 
preheating was at a = 3 in the mean medium (<5 = 0), the 
required temperature would be T ~ 7.5 x 10® K. If, on the 
other hand, heating is in the central region of a present-day 
cluster, where <5 ~ 10"^, the required temperature would be 
T ~ 2 X lO^'K. 

The time-scale of radiative cooling, defined as tcooi = 
T/T, can be written as 


tcool 


SkTfii 

2/rneA(r) 

1.0 X lO^AIagTe 


0.02 J 




(5) 


where Te = r/10®K, fi « 0.6 and /ii « 1.2 are the 
mean molecular weights per particle and per ion, respec¬ 
tively, and A _23 is the cooling function A{T) in units of 
10“^® cm®ergs“® as defined in Sutherland & Dopita (1993). 
Notice that /ti enters because Sutherland & Dopita (1993) 
defines the cooling rate per unit volume as n^niA{T), where 
Ui is the total ion density. If we approximate the cooling 
function as a power law of gas temperature, A(r) oc T~°‘, 
the cooling time in an approximately isentropic gas scales 
as 

tcooioc^--ocT“-®/Vn^(“-i/^)/®. (6) 

Thus, depending on whether a > 1/2 or a < 1/2, the 
cooling time is longer or shorter for a high-density (high- 
temperature) gas. This is different from an isothermal gas, 
where tcooi is always shorter for higher rie. 

The cooling time given above should be compared to 
the Hubble time 

tH = 3.2 X lO^ho,® (7) 


where 'H{z) = Note that 7d(z) ~ 1 

for a ^ 1 in a flat universe with Qo 0.3. As one can see, 
cooling is ineffective at low overdensities and at low pre¬ 
heating redshift {z ^ 3). In this case, the specific entropy 
S is approximately conserved (except in shocks), and the 
evolution of the temperature is given by equation (^. The 
properties of the preheated gas are then specified by its den¬ 
sity [which is proportional to (1 -I- (5)(1 + z)®] and the value 
of 5. 

One way to specify the initial entropy is to use the ob¬ 
served entropy excess in clusters and groups (Ponman, Can¬ 
non & Navarro 1999; Lloyd-Davies et al. 2000), which gives 
S ~ lOOkeVcm^. We use this as the fiducial value for 5. 
It must be emphasized, however, that in reality the value 
of S may change from place to place, but we ignore this 
complication here. 



Figure 1. Density profiles of preheated gas (assumed to be adi¬ 
abatic, with cSioo = 1) as a function of radius around dark 
haloes. Six curves are shown for redshift from 0 to 2.5 with a 
step size of 0.5. Results are shown for (a) a normal galaxy halo 
with Tvir = 10® K and (b) a dwarf galaxy halo with T^ii. = 10® K. 


In Section ^ we will discuss in some detail how the IGM 
may have been preheated. Based on the discussion there, we 
advocate a preheating process which can be summarized as 
follows; 

• The preheating was due to the starbursts (and associ¬ 
ated AGN activity) responsible for the formation of old stars 
in ellipticals and bulges. 

• The preheating was likely to be confined to protocluster 
and protogroup regions where starbursts were common. 

• Although preheating must be a continuous process and 
different regions may have been preheated at different times, 
the bulk of preheating was likely to be around a redshift 
2^2-3, where the star-formation rate and AGN number 
peak. 

The details of the preheating process are still poorly 
understood, but much of the following discussion is inde¬ 
pendent of these details. 


2.2 Accretion of the hot IGM 


Gonsider a dark halo with mass Mvir, circular velocity 14, 
virial temperature Tvir, and virial radius rvir. According to 
the spherical collapse model, these quantities are related as 
follows: 


Mvir = 




14 


10GH{z) 


WH{z) ’ 


Tvir = 


n^pVc 

2k 


( 8 ) 


The virial temperature is the same as the isothermal tem¬ 
perature for a singular isothermal sphere, but its meaning is 
less clear for more complex halo profiles. It should only be 
regarded as an approximate measure of the gas temperature 
in the halo. We model the dark matter halo as a singular 
isothermal sphere with density profile 


p(r) = 




47rGr2 


(9) 
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Figure 2. Temperature profiles of preheated gas (assumed to 
be adiabatic, with cSioo = 1) as a function of radius around dark 
haloes. Results are shown for (a) a normal galaxy halo with Tvir = 
10® K and (b) a dwarf galaxy halo with Tvir = 10® K. 


In a preheated medium where gas cooling is unimportant, 
the halo can accrete the preheated gas though adiabatic 
compression. If the gas accretion is subsonic (Tozzi, Scharf 
& Norman 2000), then we may solve for the distribution of 
gas within the halo as follows. Assuming spherical symmetry 
and hydrostatic equilibrium (the latter assumption is valid 
in a region where the sound-crossing time is much shorter 
than the age of the universe), we have 


1 dP _ GM 
p dr J .2 ’ 


( 10 ) 


where p and P are the density and pressure of the gas at 
radius r, M is the mass within radius r. For the gas outside 
the halo, we may take M to be the total mass of the halo. If 
the gas is not significantly shocked when establishing the hy¬ 
drostatic equilibrium, as is the case if the temperature of the 
preheated gas is not much lower than the virial temperature 
of the halo, the specific entropy T/rie^^ will remain at the 
preheated value S. In this case, the hydrostatic equilibrium 
equation can be solved for the gas density: 


2/3 _ 2 GMvirpnip 1 

5 kS r 


+ constant. 


( 11 ) 


Assuming that n is equal to the mean density at r —> cx3, 
and using the relation between halo mass and halo virial 
temperature, we have 


n^/3(r)=n2/® + l%!:iii, (r > r.n). ( 12 ) 

0 O T 

We caution that equation ([I^ is only approximate at large 
radii as the gas is likely not in hydrostatic equilibrium. 
Within the virial radius of an isothermal sphere, the hy¬ 
drostatic equilibrium of adiabatic gas leads to the following 
profile: 


-,2/3 (^) ^ „2/3 ("l - In — , (r < Tvir) . 

5 O \ ^vir / 


(13) 


The corresponding temperature profile is given by 


T{r) = Sn^^{r). (14) 

Some examples of the density and temperature profiles are 
shown in Figures ^ and As one can see, the density profile 
is very shallow. In fact this is the shallowest equilibrium 
profile one can get in a singular isothermal sphere, because 
any shallower profile is convectively unstable. 


2.3 X-ray emission from galactic haloes 

Given the density and temperature profiles, we can estimate 
the radiation luminosity of the hot gas in a dark halo. Since 
the temperature of the gas is about 0.1 keV for a galaxy 
halo, the radiation is in the soft X-ray band, and we can 
hope to detect such emission through X-ray observations. 
The bolometric X-ray surface brightness is given by 

‘S'x(-R) = ^ J A[T{r)]nenidy, (15) 

where r = {R^ + and R is the radius in projec¬ 

tion. The luminosity within a radius, I/x(< R), can then be 
obtained by integrating Sx over area and solid angle. Fig¬ 
ure ^ shows Sx(f?) and I/x(< R) for a present-day galactic 
halo with virial temperature Tvir = 10® K, corresponding to 
14 ~ 170kms“^. In the calculation, we have taken a metal- 
licity Z = O.IZq and neglected the contribution of the gas 
outside the virial radius of the halo as the gas tempera¬ 
ture becomes too low to be in the X-ray band (see Fig. 2). 
For simplicity, we have also assumed that all the gas is in 
the hot phase. As one can see, the predicted X-ray emis¬ 
sion depends significantly on the assumed specific entropy: 
a higher value of <Sioo leads to a less concentrated gas dis¬ 
tribution and so to lower X-ray surface brightness and lumi¬ 
nosity. For iSioo = 1, the surface brightness at R — 100 and 
20kpc are 5x ~ 2.5 x 10“® and 5 x 10“® ergs“^ cm”^sr“^, 
respectively. The total luminosity within the virial radius is 
Lx ~ 3 X 10^®ergs“^. These predicted surface brightness 
and luminosity are about a factor of 10 smaller than what 
one obtains based on the assumption that gas in galactic 
haloes is shock-heated to the virial temperature and in hy¬ 
drostatic equilibrium in the dark-halo potential (e.g. White 
& Frenk 1991; Benson et al. 2000), because of the shallower 
gas density profiles. 

It should be pointed out that the results obtained here 
are valid only for isolated late-type galaxies. For galax¬ 
ies in clusters and groups, the X-ray emission may be af¬ 
fected by the intracluster (intragroup) medium. For ellipti¬ 
cal galaxies which contain many stars, the interaction be¬ 
tween supernova-driven winds and gas accretion may also 
play an important role in the formation of X-ray haloes (e.g. 
Brighenti & Mathews 1999). Thus, the most relevant obser¬ 
vational results to compare our model predictions with are 
those obtained by Benson et al. (2000) for the diffuse X-ray 
emission from 3 late-type galaxy haloes. The upper limits on 
the bolometric luminosities of these three haloes are of the 
order 10"^^ ergs“^. These observational results are consistent 
with our model prediction, if the circular velocities of these 
haloes (at virial radius) are 14 ^ 200kms“^. If the circular 
velocities of these galaxies were as high as the measured ro¬ 
tation velocities of the disks (~ 300kms“^, corresponding 
to a virial temperature of about 0.3 keV), shock heating of 
the gas by gravitational accretion becomes significant. Cal- 
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Figure 3. The X-ray surface brightness S^{R) and luminosity 
^x(< R) (within radius R) of the hot gas around a present-day 
galactic halo with virial temperature Tvir = 10® K. Results are 
shown for three different values of Sioo- 


culations by Tozzi & Norman (2001) including this effect 
give an X-ray luminosity of ^ lO'^^ergs”^ for haloes with 
Tvir = 0.3 keV, assuming a preheated entropy similar to that 
which we adopt here. Thus, preheating can help to alleviate 
the X-ray halo problem for late-type galaxies. 


2.4 Accretion of cooled clouds 

Given the density and temperature profiles, we can esti¬ 
mate the cooling time scale once the metallicity of the gas 
is known. Figure ^ shows the ratio between the cooling time 
tcooi and the Hubble time fn assuming that the gas has a 
metallicity Z = O.IZq (the cooling time is reduced by a 
factor of about 2 if we assume Z = 0.3Zq). We used the 
cooling function given in Sutherland & Dopita (1993). No¬ 
tice that the smaller the ratio fcooi/fu, the more effective the 
cooling is. Effective cooling is expected when fcooi/fn ~ 1. 
As one can see from the figure, cooling is more effective at 
low redshift. The main reason for this is that the gas den¬ 
sity around dark haloes does not change very rapidly with 
redshift (and so the cooling time does not change rapidly) 
while the Hubble time increases with decreasing redshift. 
The situation is quite different in the absence of preheating, 
where the characteristic gas density changes with redshift as 
(1 -I- z)^ and cooling is always more effective at higher red¬ 
shift. Note that once cooling becomes important, the gas can 
no longer be considered adiabatic. Thus, the adiabatic as¬ 
sumption may break down at low redshift. The ratio Icooi/Ih 
does not change significantly with radius at a given redshift, 
implying that the cooling efficiency is quite independent of 
radius. This is very different from what one obtains based 
on the assumption that gaseous haloes are (approximately) 
isothermal spheres, where gas cooling is always inside-out. 

Because the cooling time is quite long for the preheated 
medium, most of the gas can remain in the hot phase. But 
some gas can cool due to radiative cooling. Since the pre¬ 
heated gas is expected to be clumpy, cooling is expected to 
start from regions where the gas density is enhanced, and 
a two-phase medium may develop due to thermal instabil- 



r/r - 

'' VII 


Figure 4. The ratio between cooling time and Hubble time of 
preheated gas around dark haloes. Results are shown for (a) a 
normal galaxy halo with Tvir = 10^ K and (b) a dwarf galaxy 
halo with Tvir = 10^ K. 



z 


Figure 5. The mass fraction in cooled gas (defined as the ratio 
between the mass of cooled gas Mcool Mb q, with Mb,o equal 
to /b times the total halo mass at the present time) as a function 
of redshift in protogalaxy regions which were preheated at z = 3. 
The specific entropy is assumed to be <Sioo = 1-0 the gas 
metallicity is assumed to be O.IZq. Results are shown for haloes 
with Tvir = 10®, 10® and 10® ® K. 


ity, with the hot phase at the temperature set by preheating 
and a cold phase at ^ 10“^ K (below which radiative cooling 
may be inefficient, particularly in the presence of photoion¬ 
ization heating). The cold clouds produced in this way may 
be gravitationally accreted into dark haloes. 

Since both the density and temperature of the gas may 
change with time as the dark halo grows, and since the 
medium is likely inhomogenous and multi-phase, an accu¬ 
rate estimate of the amount of cold gas that can be accreted 
in a protogalaxy can be obtained only with the use of nu¬ 
merical simulations. Here we give an approximate estimate 
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based on simple assumptions. We assnme that as a dark halo 
grows with time, its circnlar velocity K remains constant. 
For simplicity, we assume the dark halo mass profile is de¬ 
scribed by a singular isothermal sphere, and so its radius 
increases with time as rvu oc 1/H{z) (see eq. At any 
given time, most of the gas is assumed to be hot, with the 
gas density profile given by equations ® and (|^), while 
part of the gas in high-density regions cools. Assuming fur¬ 
ther that the total gas mass within a spherical mass shell 
does not change with time, we can estimate (at any given 
redshift) the cooling time using the equilibrium density and 
temperature of the mass shell in consideration. From the 
definition of the cooling time, the fraction of the gas that 
has cooled by redshift z can be estimated from 

/cool(z) = f 

J tcool 

X f ^^{l + 5)[l-Uoi{z')]{l + z'fdz', (16) 

•f 2ph ° 

where Zph is the redshift at which the gas is preheated, and 
the factor [1 — fcooi{z')] in the integrand takes into account 
the depletion of the hot gas due to cooling. Since the density 
of a fixed mass shell is known as a function of a in our 
model, we can obtain fcooi{z) for each mass shell. We can 
then integrate over all the mass shells to obtain the total 
amount of gas that has cooled in the protogalaxy region by 
some redshift 2 . Figure]^ shows the ratio, Mcooi/Afe.o (where 
Mb,o is equal to /b times the total halo mass at the present 
time) as a function of 2 . Several important conclusions can 
be drawn from the figure, (i) In all cases, significant cooling 
occurs at low redshift, for the reasons discussed earlier in this 
subsection, (ii) For small haloes, rapid gas cooling occurs at 
2^1 when the temperature of the preheated gas reaches 
the level T ^ 2 x 10® K where the cooling rate peaks. In this 
case, the large amount of cooled gas may trigger episodes 
of star formation, which may in turn drive outflows and 
prevent further accretion of cooled gas. (hi) For large haloes, 
the fraction Mcooi/Mb,o [which is roughly proportional to 
(1 -I- 2 )“®/^ at low 2 ] is signihcantly smaller than 1, and so 
only a fraction of the gas in the protogalaxy region can cool 
and be accreted. 


3 GALAXY FORMATION IN A PREHEATED 

INTERGALACTIC MEDIUM 

In the last section we have shown that gas accretion by 
galaxy haloes in a preheated medium has several distinct 
properties. In this section, we discuss the implications of 
such gas accretion for the formation of galaxies. 

Based on the discussion in the last section, we may 
roughly divide stars in a galaxy into two populations: one 
formed in starbursts before preheating and the other formed 
due to the accretion of gas in the preheated medium. The 
relative importance of these two populations may determine 
the morphological type of the galaxy, because the early star- 
bursts were likely to form a bulge component, while an ex¬ 
tended disk can only form through gas accretion. 


3.1 The formation of disk galaxies 

As we have seen in Section the gas accretion in a preheated 
medium by dark haloes occurs at rather low redshifts (more 
than half of the gas is accreted at 2 1). Thus, galaxy disks 

are expected to form late in this model. In fact, late forma¬ 
tion of present-day galaxy disks is required in the standard 
model of disk formation (where disk size is determined by 
the specific angular momentum of dark haloes), because the 
formation at high redshifts (> 1) predicts disks that are too 
compact to match the observed sizes (e.g. MMW). The late 
formation of galaxy disks (relative to the bulges) may also 
be required in order to solve the G-dwarf problem (Ostriker 
& Thuan 1975). 

Another important property of the disk formation in a 
preheated medium is that only part of the gas in a proto¬ 
galaxy region can settle into a disk in the galaxy halo. Based 
on the discussion in Section |^, this fraction is about / = 0.3 
- 0.5 for galaxy haloes, and so the ratio between disk mass 
and halo mass is ^ /flB.o/Ho ~ 0.04 - 0.06 (where 
we have assumed Hb.o = 0.02h~^, flo = 0.3 and h — 0.7). 
This fraction is in fact required in the standard model of 
disk formation. Indeed, if all of the gas in a protogalaxy 
region can settle in the disk, then md ~ 0.13 in the standard 
ACDM model. Such a high value of irid leads to too massive 
disks which have too peaked rotation curves and may be 
unstable (MMW). The model present here naturally leads 
to the required reduction. 

In the standard model of disk formation, the specific 
angular momentum of disk material is assumed to be com¬ 
parable to that of dark matter. Since disk material and dark 
matter is well mixed in a protogalaxy, this assumption may 
be valid before the collapse of the dark halo. However, in or¬ 
der for the assumption to hold for the formed disk, one has 
to make the assumption that there is no signihcant angular- 
momentum transfer from the disk material to dark matter 
during the formation of disks. As shown by the simulations 
of Weil et al. (1998), this requires the gas to remain in a 
diffuse form before it collapses into the final halo. This is 
exactly the case in the present model, where most of the 
gas remains in a hot, diffuse form until late time (2 ^ 1) 
when the final halo collapses. Since this gas does not suffer 
the same dynamical friction as dark matter clumps, the final 
distribution of specific angular momentum of the gas may 
be different from that of the dark matter. 

To see this point more clearly, we consider a simple 
case where a smaller halo, Mvir, merges into a larger halo. 
For simplicity, we assume the total mass profiles in 
both haloes are described by that of a singular isothermal 
sphere, while the gas distribution in the small halo is given 
by eq. (p^. The smaller halo is taken to be always on a 
circular orbit as it sinks towards the center of the larger 
one due to dynamic friction. Because of tidal stripping, the 
outer part of the smaller halo will be truncated, and we 
estimate the truncating radius rt on a orbit of radius r' by 
p(< rt) = p'{< r'). For the assumed density prohle, we have 
rt = (»'vir/r(,;,.)r', and the remaining mass of the smaller halo 
within rt is M(< rt) = V^rtfG. For material on circular 
orbits in the bigger halo, the specihc angular momentum is 
j(r') ~ V^r'. It is then easy to show that after merger the 
mass of dark matter that originated in Mvir with specific 
angular momentum smaller than j is 
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Figure 6. The fraction of the total halo mass M in progenitors with circular velocities exceeding Vc as a function of redshift. Panel (a) 
shows the results for a present-day rich cluster with mass M = 10^^ Mq, while panel (b) shows the results for a present-day group with 
mass M = 2.5 x The circular velocities of progenitors are labeled on the curves in units of kms~^. The calculations are for the 

standard ACDM model based on the extended Press-Schechter formalism (e.g. Bond et al. 1991). 


(17) 

>ir 

where jvir = For the gas, the corresponding mass is 

given by 

ri 

Mgas(<j)=/ 47rpgas(r)r^dr, (18) 

Jo 


where Pgas(r) is the gas mass density around the smaller 
halo, and Vj = rviij/jvii- For the gas density profile given 
by equation di), 


A7gas(< j) 


47rPgas>'vir(l + «) 




a 

1 + a 


3/2 


In* 


d*. 


(19) 


where a = (4/5)rvir/(5n^^^). Since the logarithmic term 
changes only slowly with x, we have, to a good approxima¬ 
tion, 


lUgas (< i) OC 



( 20 ) 


Numerically, we found that for a = 0.1, Mgas(< j) oc 
and for a — 10, Mgas(< §) oc . Clearly this distribution 
is much steeper than that given by equation (§) for the 
dark matter. In fact, the shallower the gas density profile 
is, the steeper the dependence of Mgas(< j) on j. Phys¬ 
ically, because of the shallow density profile, most of the 
gas is stripped in the outer part of the bigger halo and so 
retains higher specific angular momentum; only a small frac¬ 
tion of the gas can sink to the inner part and lose a large 
amount of angular momentum. This effect exists wherever 
the gas distribution is more extended than the dark matter. 


We emphasize that this effect remains even if we have mul¬ 
tiple accretion events that occur with different angles and 
this will lead to a reduction in the net specific angular mo¬ 
mentum of the accreted gas. However, the specific angular 
momentum of the gas will still be higher than that of the 
dark matter, if both are reduced by the same factor. This 
may help to solve the problem pointed out by Bullock et al. 
(2001) and van den Bosch et al. (2001) that the assumption 
that the gas and dark matter in dark haloes have the same 
angular-momentum distribution leads to too concentrated 
disks. 


3.2 Formation of cluster galaxies and 
morphological segregations 

As one can see from Figure |^, significant fraction of the to¬ 
tal mass of a present-day rich cluster is already in group¬ 
sized systems with circular velocity 14 ^ 400kms“^ at 
z ^ 2. These groups are expected to be associated with high- 
density peaks in a protocluster region, and so are likely to 
end up in the central part of a cluster. Since galaxy haloes 
cannot have accreted significant amount of gas from the pre¬ 
heated gas at 2 > 1 and since cooling of gas is already inef¬ 
ficient in such group haloes, we expect that galaxies which 
end up in the central part of clusters contain mainly old stars 
that have formed through starbursts before preheating. As 
the haloes of these galaxies merge, the stellar components 
also merge due to dynamical friction to form bigger galax¬ 
ies. Such mergers of galaxies continue until the galaxies are 
incorporated into clusters in which dynamical friction is no 
longer effective. We identify this to be the formation path 
for giant ellipticals in clusters. This formation is in fact sim- 
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ilar to that in the standard model where cluster ellipticals 
are formed by mergers of galaxies (e.g. Kauffmann 1996). 
The effect of preheating is to suppress the formation of new 
stars, because gas accretion was prohibited after preheating. 
Since a group halo at 2 ; ^ 1 contains many smaller haloes 
at the time of preheating, the formation of a giant ellipti¬ 
cal is expected to be a result of multiple mergers of stellar 
systems formed through early starbursts, and clusters ellipti¬ 
cals are therefore expected to contain mainly old stars. This 
early formation of elliptical galaxies (at least their stars) in 
gas-rich systems may be responsible for their tight colour- 
magnitude relation (e.g. Bower et al. 1992) and other proper¬ 
ties which are distinct from present-day spiral galaxies (e.g. 
Ostriker 1980). 

As one can also see from Figure ^ about 30% of the to¬ 
tal mass of a present-day rich cluster remains in small haloes 
{Vc ^ 200kms“^) at z ^ 1. This fraction is expected to end 
up in the outer part of the cluster. These haloes can accrete 
significant amount of gas before they are incorporated into 
larger systems. The accreted gas may form new stars on a 
disk around an old bulge, although young stars may also be 
added to the bulge either due to disk instability or due to the 
accretion of satellites containing young stellar component. If 
such systems merge into groups before being incorporated 
into the final cluster, the galaxies they contain may also 
merge, producing ellipticals with young stellar populations; 
otherwise they remain as spiral or SO galaxies if their disks 
are not destroyed later in the cluster environment. 

Similar discussion can be made for galaxy formation in 
smaller systems, such as present-day groups of galaxies. The 
only difference here is that, by definition, galaxy haloes can 
survive to later times before being incorporated into group 
haloes with K ^ 300kms“^ (see Fig. |b) and so can accrete 
more gas to form disks. The mergers of such disks may also 
form elliptical galaxies, but unlike giant ellipticals in clus¬ 
ters, these ellipticals should contain a lot of younger stars. 
Since the merger progenitors here contain gaseous disks, 
such mergers can produce the tidal tails and starbursts ob¬ 
served at low redshifts. 

If we use the disk-to-bulge ratio to represent the mor¬ 
phological type of a galaxy, the above discussion suggests 
that the gas content, stellar age and local environment of 
galaxies should change systematically along the morphology 
sequence. Late-type galaxies are formed in low-density re¬ 
gions where gas accretion lasts for a longer period, so they 
are gas richer, bluer and residing in lower density than early- 
type galaxies. This trend is qualitatively consistent with ob¬ 
servations (e.g. Roberts & Haynes 1994; Dressier 1980; Post¬ 
man & Geller 1984). However, quantitatively this prediction 
has to be studied taking into account of dynamical evolution 
and merging between galaxies. Such a detailed calculation 
in a cosmological context is clearly needed but is beyond the 
scope of this paper. 


4 PREHEATING THE INTERGALACTIC 
MEDIUM 

The preheating of the IGM has been discussed quite ex¬ 
tensively in connection to the observed entropy excess in 
clusters and groups (Kaiser 1991; Evrard & Henry 1991; 
Ponman et al. 1999; Tozzi & Norman 2001 and references 



Figure 7. The mass fraction of the universe in haloes with circu¬ 
lar velocities exceeding Vc (labeled on curves in units of kms“^) 
as a function of redshift. The calculations use the halo mass func¬ 
tion given in Sheth, Mo & Tormen (2001), and are for the stan¬ 
dard ACDM cosmogony with SIq = 0.3, 12^ = 0.7, h = 0.7 and 
as = 0.9. The thick dashed curve shows the mass fraction in 
haloes which can trap photoionized gas [adopted from Gnedin 
( 2000 )]. 

therein). In this section we discuss further some points which 
are closely related to our discussion of galaxy formation. 

4.1 Preheating by starbursts and AGNs 

If we assume a stellar mass function (IMF) of the standard 
Salpeter form [n{M)dM oc for O.IMq < M < 

50M©] and that each star of mass greater than 8 Mq releases 
lO^^Esierg in kinetic energy in a supernova explosion, the 
total energy output (Esn) is related to the total amount of 
stars formed (M*) by 

« 7.5 X lO'‘®E5i(M*/M0)erg. (21) 

If a fraction e of this energy is to heat gas with total mass 
M, each gas particle will gain an energy (in terms of tem¬ 
perature) : 

Tsn - 1.8 X lO^eEsi K. (22) 

The ratio between this temperature and that given by equa¬ 
tion (|) is 

^ ^ 4 IOCE 51 5ro'o(l + 5)-"^"(l + 2 )-" . (23) 

If we take 5 = 0, 2 = 3 and <Sioo = 1, then eM*/M ~ 0.04. 

X-ray observations of galaxy clusters show that the 
mass in stars is about 10% of the total mass in the X-ray 
gas (e.g. Bohringer 1995; Balogh et al. 2001). In this case 
M^,/M ~ 0.1. A similar number can be obtained as follows. 
The observed mass density in stars is fl* ~ 2 x 10“®h“'^ 
and one-third to a half of this may be in old stellar sys¬ 
tems, such as ellipticals and bulges (e.g. Fukugita, Hogan 
& Peebles 1998). If we take the standard ACDM model. 
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the fraction of cosmic mass in protogalaxy regions (i.e. re¬ 
gions which form haloes with circular speeds ~ 200kms“^ 
at the present time) is about 0.3 (see Figure ^). Thus 
Af*,oid/A/ ~ fl*,oid/(0.3flB,o) ~ 0.1. If the heating of the 
IGM were all due to star formation, an efficiency e ~ 0.4 
would be required. This efficiency is higher than usually as¬ 
sumed (e.g. Valageas & Silk 1999; Wu, Fabian & Nulsen 
2000; Kravtsov & Yepes 2000). However, an efficiency as 
high as this is in fact indicated by observations of supernova- 
driven winds in local starburst galaxies (Heckman et al. 
2000). Thus, if most stars at high redshifts were formed in 
systems reminiscent of local starbursts, as suggested by the 
observed properties (high star-formation rate, compact size) 
of high-redshift galaxies (see Heckman 2001 and references 
therein), the stellar energy sources may be sufficient to pre¬ 
heat the IGM to the required level. More quantitatively, the 
observed mass-outflow rate in a local starburst is usually 
comparable to the star-formation rate, and the wind veloc¬ 
ity is typically 500 - 600 kms“^. If the kinetic energy of the 
wind is to be thermalized in the IGM, it can heat as much 
as ten times Mwind (the mass of the wind) of the IGM to 
the required temperature at 2 ~ 3 [see equation (^)]. Since 
Mwind ~ M*, the total preheated mass is M ~ lOM*, con¬ 
sistent with the star/gas mass ratio estimated above. 

An alternative energy source for the preheating may 
be AGNs. For example, radio galaxies can provide a total 
energy output about two orders of magnitude larger than 
stellar sources, and so have the potential to heat the IGM 
to the required level (e.g. Inoue & Sasaki 2001). Since the 
mass of central black hole (which powers the AGN) is found 
to be proportional to the bulge mass in a galaxy (e.g. Magor- 
rian et al. 1998; Gebhardt et al. 2000), the energy output 
from AGNs in a region large enough to contain many AGNs 
should be roughly proportional to the energy output from 
stellar sources. Thus, the inclusion of additional energy out¬ 
put from AGNs is equivalent to increasing the value of e. 
The uncertainty here is again how large a fraction of the 
total energy released by AGNs can be thermalized in the 
IGM. 

4.2 The epoch of preheating 

If the preheating was indeed done by activity associated with 
star formation and AGNs, the likely epoch for preheating 
is when the star-formation and AGN activity peaks. The 
observed comoving number density of AGNs peaks around 
2 = 2 to 3 (Shaver et al. 1996). The star formation history in 
the universe can be inferred from the observed H*( 2 ) which 
is the change rate of the density parameter of stars. Ob¬ 
servations showed that Cl^, remains approximately constant 
(~ 10“® Gyr“^ ) in the redshift range 2 = 2-4, and de¬ 
creases rapidly with decreasing 2 at 2 < 1 (e.g. Blain et al. 
1999). For the ACDM model, the age of the universe at 2 ~ 3 
is about 1.5 Gyr, and so the mass density parameter of all 
stars formed before 2 = 3 is ~ 10“®, consistent with the 
total amount of stars observed in old stellar systems. Thus, 
if the mass of gas that can be heated up to the required 
entropy is about ten times that of the old stars, the den¬ 
sity parameter of the heated gas is Hgas ~ 0.01 ~ 0.2flB,o, 
which is similar to the gas mass in proto-group regions (see 
Fig. 0 ) . From these considerations, we see that the epoch of 
preheating is likely to be around 2 = 2 to 3. 


Similar conclusion can be drawn from theoretical con¬ 
siderations. As shown in Figure about 20% of the total 
mass of a present-day cluster is already in progenitors with 
circular velocities 14 ^ 200kms“^ at 2 ~ 3. This fraction is 
only slightly smaller for a present-day group (see Fig. |b). 
Most systems that have formed by redshift 3 will suffer a 
major merger in a Hubble time (see Fig.14 in Shu, Mao 
& Mo 2001), and so most of them may form stars in star- 
bursts. If half of the gas in these systems forms stars, and 
the other half escapes as outflows with the escaping velocity, 
Uesc ~ 314 [where we assumed a singular isothermal sphere 
and the gas escapes at roughly one disk scale length, cf. 
equation (12) in MMW and equation (2-192) in Binney & 
Tremaine (1987)[, then the IGM in the protogroup (or pro¬ 
tocluster) regions can be preheated to the required level at 
2 ~ 3. As one can see from Fig |^, about half of the total 
mass of a present-day cluster is in progenitors with circular 
velocities 14 ^ 50kms“^ at 2 ~ 3. Since such haloes cannot 
trap much gas in the presence of the general UV background 
(e.g. Efstathiou 1992; Gnedin 2000), more than half of the 
protocluster (protogroup) gas must be in diffuse form. Thus, 
most of the gas to be heated is in a diffuse medium outside 
collapsed haloes. The situation here is different from that 
in previous feedback schemes without preheating, where the 
goal is to heat up and get rid of the gas that has already 
been assembled into galaxies or galactic haloes. 

In reality, the time of preheating must be different in 
different regions. In high-density regions where a significant 
fraction of the mass was already in galactic haloes at high 
redshifts, the time of preheating is expected to be early, 
while in low-density regions where galaxy haloes form at 
low redshifts, the time of preheating is expected to be late. 
In regions where starbursts never occur, the IGM is not pre¬ 
heated. A detailed modelling of how and when different re¬ 
gions in the universe were heated up depends on how the 
heating sources form in the cosmic density field and how 
their energy is thermalized in the IGM. This is a compli¬ 
cated problem and will not be considered in this paper. 

4.3 Preserving the Lyo forest 

If all of the IGM were heated to the temperature implied 
above, then there would be no photoionized gas to produce 
the observed Lyo forest. This is not allowed, because the 
Lyo forest is observed. However, the above discussion does 
not require the preheating to be everywhere in the IGM. In 
fact, in order to explain the observed entropy excess, it is 
only necessary to preheat the gas which ends up in present- 
day groups and clusters. Since groups and clusters are places 
where vigorous star-formation and AGN activities are ex¬ 
pected in the early time, it is also likely that preheating oc¬ 
curred mainly in protocluster (protogroup) regions. In fact, 
star-forming galaxies at high redshifts, such as Lyman-break 
galaxies and sub-mm sources observed at 2 ~ 3, are expected 
to be located in high density regions which will eventually be 
in present-day groups and clusters (e.g. Mo, Mao & White 
1999, and references therein). It is therefore possible that the 
IGM in the low density regions, which may be responsible for 
the observed Lya forest, is not affected significantly, while 
the high density regions are heated. To demonstrate that 
this is likely the case, we can estimate the fraction of the to¬ 
tal mass which is in systems more massive than present-day 
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groups (i.e. with circular speeds ^ 400kms“^). As shown 
in Figure ^ this fraction is about 20% in standard ACDM. 
In fact more than half of the total mass is in systems which 
are too small to trap photoionized gas. These systems are 
outside protogalaxy regions, and so may not be influenced 
significantly if preheating is confined to proto-group regions 
by the associated gravitational potentials. 

Of course, preheating may produce some effects that 
can be observed in the Lya forest (Theuns, Mo & Shaye 
2001). If the preheating is due to supernova-driven winds 
from starbursts, some of the metals ejected by the heating 
sources may leak to the Lyman forest and contaminate it, 
although most of the metals may be hidden in the hot phase. 
This may in fact be the origin of the ‘missing metal’ problem 
that the total amount of metals observed at 2 ~ 3 seems 
smaller than that expected from the observed star formation 
history (e.g. Pettini 1999). 


5 DISCUSSION AND SUMMARY 

In this paper, we have outlined a scenario of star-formation 
feedback in which the IGM in galaxy-forming regions was 
preheated. The preheating of the intergalactic medium in 
high-density regions of the cosmic density field may be re¬ 
sponsible for the entropy excess observed in low-mass clus¬ 
ters of galaxies and yet may not destroy the bulk of the Lyo? 
forest. The subsequent galaxy formation in preheated media 
is different from that in models without preheating, because 
the gas component no longer follows the bottom-up clus¬ 
tering hierarchy of the dark matter component on galaxy 
scales. Because of the initial entropy, the hot gas around a 
galaxy halo has a very shallow profile, and so emits only 
weakly in X-ray. This may be the reason why extended X- 
ray haloes are not yet detected around spiral galaxies. Ra¬ 
diative cooling of the preheated gas around galaxies may 
become more significant at lower redshift when the cosmic 
time becomes longer. Unlike in the absence of preheating, 
where cooling is always inside-out in a halo, the cooling ef¬ 
ficiency in a preheated halo does not change rapidly with 
radius. Because of thermal instability, the cooled gas may 
form cold clouds which can then sink toward the centre in 
the potential wells of galaxy haloes to form galaxy disks. 
The total amount of the gas that a galaxy halo can accrete 
by the present time is significantly smaller than the total 
gas in the protogalaxy region. Since most gas is accreted in 
diffuse clouds, the gas may not lose angular momentum to 
dark matter due to dynamical friction. Cluster ellipticals are 
produced by the mergers of stellar systems formed prior to 
the preheating, because these systems cannot accrete much 
gas from the preheated medium to form new stars before 
they are incorporated into large haloes where gas cooling is 
inefficient. Large galaxy disks form in low-density environ¬ 
ments where accretion of gas can continue to the present 
time. Mergers of such disks may trigger new starbursts and 
form field ellipticals. Thus, the stellar age, gas content, and 
local environment of galaxies are expected to change sys¬ 
tematically along the morphology sequence. 

The scenario proposed here can be tested by observa¬ 
tions of its assumptions and predictions. If part of the IGM 
was indeed heated up to a high temperature, we may hope to 
find some signatures of such heating in the Lya forest (The¬ 


uns, Mo & Schaye 2001). Since the preheating is most likely 
in regions of vigorous star formation, observations of the 
Lya forest near star-forming galaxies, such as the Lyman- 
break population (e.g. Steidel et al. 1999), are important for 
probing the properties of the preheating. Because the pre¬ 
heated gas contains metals and is at a temperature much 
higher than 10"* K, observations of absorption lines of colli- 
sionally ionized metal species (e.g. O VI) are also important. 
Another test may come from observations of the diffuse X- 
ray emission from spiral-galaxy haloes. As discussed in Sec¬ 
tion 1^, current observational results are still too uncertain to 
provide stringent constraints. Future observations from the 
more sensitive detectors, such as XMM and Ghandra may 
be used to probe the gaseous haloes predicted here. 

Yet another test of the model may come from observa¬ 
tions of the cold gas around galaxy haloes. An important 
probe here is provided by QSO absorption line systems as¬ 
sociated with low-redshift galaxies. Such systems have been 
observed. For example, normal galaxies at moderate red- 
shifts (z ^ 1) are observed to produce strong Mgll systems 
whenever the impact parameter to a QSO sightline is smaller 
than ~ 40/i“^ kpc (Steidel 1995 and references therein), sug¬ 
gesting that each galaxy may possess a halo of cold clouds. 
Weaker Mgll systems are also observed to be associated with 
galaxies at larger impact parameters (e.g. Rigby, Charlton & 
Churchill 2002). These absorption systems may be produced 
by the cooling clouds predicted in the present model. Indeed, 
there are theoretical investigations along this line (e.g. Mo 
1994; Mo & Miralda-Escude 1996; Lin et al. 2000). Assuming 
that strong Mgll systems are produced by gas clouds cool¬ 
ing from and pressure-confined by hot haloes similar to those 
considered here. Mo & Miralda-Escude found that many of 
the observed properties of the absorber/galaxy systems can 
be reproduced. Clearly, more observations of this kind can 
give stringent constraints on the gaseous galaxy haloes we 
are proposing here. 

On the theoretical side, there are several important is¬ 
sues for which further investigations are required. For the 
preheating, we need a physical model to predict how the 
IGM was heated. For simplicity, we have modelled the pre¬ 
heated IGM as a uniform medium. This is clearly a simpli¬ 
fication, as the density contrast is likely a smooth function 
of radius, decreasing from the center of the preheating en¬ 
ergy sources (bulges or AGNs) to the preheated IGM. Also 
if the preheating is accompanied by outflows, the interplay 
between outflows and gas accretion needs to be modelled 
in detail. Some of these issues are best addressed using nu¬ 
merical simulations. Finally, we also need to incorporate the 
present scenario of galaxy formation into realistic cosmolog¬ 
ical models to make detailed predictions for the properties 
of the galaxy population. We intend to return to some of 
these issues in future papers. 
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